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Abstract 



It was recently observed that the one dimensional half-filled Hubbard model reproduces 
the known part of the perturbative spectrum of planar M = 4 super Yang-Mills in the 
SU(2) sector. Assuming that this identification is valid beyond perturbation theory, we 
investigate the behavior of this spectrum as the 't Hooft parameter A becomes large. We 
show that the full dimension A of the Konishi superpartner is the solution of a sixth order 
polynomial while A for a bare dimension 5 operator is the solution of a cubic. In both 
cases the equations can be solved easily as a series expansion for both small and large 
A and the equations can be inverted to express A as an explicit function of A. We then 
consider more general operators and show how A depends on A in the strong coupling 
limit. We are also able to distinguish those states in the Hubbard model which correspond 
to the gauge invariant operators for all values of A. Finally, we compare our results with 
known results for strings on AdS^ x S 5 , where we find agreement for a range of i?-charges. 



1 On leave from Institutioncn for Teoretisk Fysik, Uppsala Universitet, Uppsala, Sweden 



1 



1 Introduction 



The underlying integrability of planar M = 4 super Yang-Mills is expected to be crucial 
in better understanding the AdS/CFT correspondence [U 121 El- At one loop integrability 
is well-established, since the gauge invariant operators can be mapped to the states of 
an integrable spin chain jU In this map the dilatation operator plays the role of 
the spin chain Hamiltonian. But higher loop integrability requires the inclusion of all 
loops, making it difficult to prove. Nevertheless, important progress has been made in 
establishing higher loop integrability, starting with the work in j^j, and many further 
tests have been carried out |H] • The dilatation operator is known with some amount of 
confidence up to five loops in the reduced 577(2 1 3) sector but an all loop operator seems 
out of reach, or if it could be reached, not particularly useful. 

On the other hand, it might not be necessary to find the dilatation operator if one 
knows the 2-body S-matrix and a dispersion relation for the "magnons" on the spin 
chain. With the 5-matrix, one can construct Bethe ansatz equations for the magnon 
momenta and from there compute the energies of the states QDj. In jllj . Beisert 
Dipple and Staudaucher (BDS) gave a conjecture for the asymptotic Bethe ansatz in the 
SU(2) sector which is consistent with the three loop dilatation operator and which has 
a form similar to an integrable equation derived from classical string motion on R x S 3 
[T2] . Unlike a previous conjecture [T5j . the BDS Bethe equations are consistent with all 
loop BMN scaling Using the BDS equations it is possible to compute the largest 

energy state on the chain [TniUS], where it was discovered by Rej, Serban and Staudacher 
(RSS) that it had precisely the same energy as a Lieb-Wu state for the one-dimensional 
half-filled Hubbard model ^7], with a hopping parameter dependent on the 't Hooft 
parameter A [13]. RSS were also able to show that integrating out the hopping terms for 
the half filled states of the Hubbard model led to an effective Hamiltonian that matched 
the 3 loop dilatation operator. They further showed that the magnon dispersion relation 
and the asymptotic S'-matrix for two magnons reproduces the BDS result. Finally, RSS 
demonstrated that the Hubbard model's Bethe equations could be used to find a series 
expansion for operators with small bare dimensions, including the dimension 4 operator 
that is part of the Konishi multiplet. 

It is not presently clear why the Hubbard model should reproduce perturbation theory. 
This is an important question that hopefully can be answered. In the meantime, one can 
treat the RSS identification of gauge invariant operators with a subset of Hubbard model 
states as an ansatz, even at strong coupling, and explore its consequences. 
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While Bethe ansatz equations are in principle solvable, actually finding the solutions 
is often cumbersome. Even the 4 and 5 site series solutions in [15J are quite involved. 
However, the Hubbard model has a simple Hamiltonian, and since it only acts on a finite 
system of states, it is possible to find the desired energies by solving a polynomial equation. 
For example, there are 36 half-filled states with an equal number of up and down spins on 
the 4 site Hubbard model. But invoking symmetries reduces the number of relevant states 
to 6. On the half-filled 5 site Hubbard model with three up and two down spins there are 
100 states, but symmetries reduce these to nine. However, it turns out that the resulting 
ninth order polynomial is itself factorizable into a cubic and a sextic polynomial, where 
the solution corresponding to the operator in the gauge theory is a zero of the cubic. 

A key point about the Hubbard model is that for large A it reduces to a theory of free 
lattice fermions. However, one still needs to determine how the weakly coupled gauge 
theory states flow to the lattice fermion states as the coupling is increased. We will show 
how to do this explicitly for a particular class of states using the nested Bethe ansatz and 
sketch the argument for more general states. Along the way we will give a criteria that 
identifies which Hubbard model states are to be identified with gauge invariant operators. 

Once we are in the strong coupling region where A >> 1, we can compare the results 
to known results in string theory. In particular, we will consider the case of operators with 
bare dimension L and i?-charge L — 2, but which are outside the BMN scaling range [14]. 
We will show that there is a range of .R-charges where the results agree with a prediction 
in [2] for the dimension of these operators. 

In section 2 we consider specific solutions for Hubbard chains with 4 and 5 sites 
respectively. We identify the states that correspond to gauge invariant operators and 
show that A is a solution to a polynomial equation. In section 3 we consider more general 
operators and compute A for these operators in the limit of large 't Hooft coupling A. In 
section 4 we compare our results from section 3 with string theory results. 



2 Solutions for 4 and 5 sites 

In this section we compute the dimensions for operators which are not chiral primaries, and 
that have bare dimensions 4 and 5, using the Hubbard model ansatz. The Hamiltonian 
for the Hubbard model which reproduces the BDS Bethe ansatz is given by ^S] 




H = TJn I E K^+i + c Wi c ^) + eH, L c.,i + e-^cl^L ) + ^(1 - 4 jiCT , iC [. Ciji ) 

1=1 

(2-1) 



3 



where (f) = ir/2 (0) if L is even (odd). For what follows we will leave the phase <fi arbitrary, 
eventually setting it to the values used in [T3]. The parameter g is given by 



2= (gy M ?N = X 
y Stt 2 8tt 2 V 1 



We have chosen to normalize the Hamiltonian as in ^8] so that the energy is given by the 
bare dimension when g = 0. The parameter <ft arises from an Aharanov-Bohm flux which 
is needed to reproduce the BDS results [T^]. The last term in (|2.1|) is an interaction term 
between the lattice fermions which lowers the energy if two fermions are occupying the 
same site. The states corresponding to the gauge invariant operators are not the lowest 
energy states The Hamiltonian is also invariant under the shift symmetry 

r ■ — > p i( t>/ L r • , -, r 1 — > p~ i< l'/ L J n <- T 

c ff , L - e^ +1 )/V,i < L - e^ 1 ^ . (2.3) 

Let us first consider the Hubbard model with just 4 sites. There are 36 possible 
half-filled states and the Hamiltonian can mix many of these states amongst each other. 
However, the symmetries of the system will reduce the mixing to a smaller subset of 
states. Since the total spin is conserved, we can reduce our space to states with total spin 
zero. We can further reduce the states to those that are invariant under the shift in ()2.3|) . 
A basis for these states is given by 

\A) = -^=(2|UU> + 2|iTiT>-|tTU>-|tUT>-|UTT>-|im>) 



2Vs 

\Bi) = ^|(||0U> + U|0T> + |U|0)-e^|0U|> 



Io|T>-|Tlo|)-||T|o) + e^|oiTl)) 



\B 2 ) = ^(e^ /2 | It0l) + e-^ 2 | |lT0>-e^ 2 |0|lT)-e^ 2 | T0|I> 
-e~^ /2 | U T > - e~^ 2 \ Til 0) + e^ 2 |0 ]\i > + ^ /2 \ 1 Tt >) 



\B 3 ) = ^=(|0lTi> + |i0|t> + |Ti0|)-e-^|:U0) (2.4) 

— |o XXT > — I T XI > - I it 1 > + e-**| tit o>) 

|Ci> = ^(e^ /2 | I0l0)-e^ 2 |0l0i>) 

|c 2 > = \(\ x 00 : > -e-^i noo> + |ono>-e^|oon>) 
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The states are defined by the fermion occupancy at each site, with the fermion creation 
operators ordered by increasing lattice number and with double occupied sites ordered 
with the down creation operator to the right of the up. So for example 



(2.5) 



I not) = c { jlC { j2 4 )4 c{ )4 |o>. 

Using the basis of states in (|2.4j) it is straightforward to show that H in (j2.1j) is given by 





H 



( 4 V3g 
V3g 3 

e~^ /2 
V3g 

e-^ 2 V2g 

\ g 



e^ 2 V2g 



< -V2g 



-e-^ 2 V2g 








V3g 


-e^ 2 V2g 
3 

-e^ 2 V2g 
9 



\ 

e^ 2 V2g g 



-e-^ 2 V2g g 

2 

2 / 



(2.6) 



Hence, the secular equation giving the eigenvalues is a sixth order polynomial and is given 
by 



det(A-F) 





A 6 - 17A 5 + (119 - 16# 2 )A 4 - (439 - 176# 2 )A 3 

+ (900 - 716g 2 + 64 cos 2 - g 4 )A 2 - (972 + 1276# 2 - 320 cos 2 | g 4 )A 



+432 - 840# 2 + 400 cos 2 | g A . 



(2.7) 



One can easily solve this as a series expansion, where the relevant solution is given by 
A 



4 + 6 # 2 -6(2 + 3cos0) # 4 + 6(7 + 45 cos 0) g 6 
-6(53 + 486 cos <p + 144 cos 2 <p) g 8 + 6(754 + 4515 cos <\> + 4815 cos 2 . 
-6(10631 + 42597 cos + 94905 cos 2 <p + 11637 cos 3 <p) g 12 
+12(65327 + 243243 cos <p + 718704 cos 2 + 299754 cos 3 0) g 14 + . . . 
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{21 



which matches the expansion found in ^Sj after setting = 7r/2. 

Since the polynomial in ()2.7|) is quadratic in g 2 , we can write a simple closed form 
solution for A in terms of A, 



A 



rr 



cos 2 0/2 



(A -2)(A-3) 
2A - 5 



2A - 7 - \\ (2A - 7) 2 sin 2 - + cos 2 - , (2.9) 
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which after setting <p = tt/2 can be expressed in terms of the dimension of the Konishi 
operator, A^- = A — 2, as 

A = 2tt 2 A ^_~ 1} (2A g - 3 - ^2A\ - 6A y + 5) ■ (2.10) 

From the expression in ()2.9|) it is easy to see that the behavior of A for large A is 

Jl + sin 0/2 r- 

A _v z-^v 7 !. (2.11) 

7T 

A similar expression with a slightly higher coefficient was found in j!9| |2*U] . where the 
dilatation operator was reduced to the linear part of the spin-spin interaction. 

We can also see the behavior in (|2.11|) arise from the original Hamiltonian in (|2.1|) . 
where in the limit of large coupling, the hopping term will dominate over the contact 
term. In this case the theory reduces to free lattice fermions with a dispersion relation 

m 

\/\ 

e(k) = j— cos(k + 4>/L) (2.12) 
2tt 

where k = 27m/L. Hence with two up fermions and two down fermions we can reproduce 
(j2.11|) by having one fermion of each type with k = and — 7r/2. 

It turns out that the unprotected dimension 5 operator has an even simpler solution. 
One starts by choosing a basis of states with total spin 1/2. The possible states that can 
mix with each other are those invariant under shifts of the lattice sites. The Hubbard 
Hamiltonian is also invariant under parity transformations that map site i to L — i. Since 
the states with all sites singly occupied are invariant under this transformation j^J, we 
only need to consider states that are parity even. This then leaves the following nine 
states that can mix with each other 

\A) = -^(| TTTII ) c - | TTITI 

\Bi) = ^(|Io|Tt> c -||om> c +|oUTT> c -|o|m> c 

\B 2 ) = J-(|0t|lt> c +|0t|U> c -2|0||tt> 



2V3 



C ' I 1 1 *♦* IG I I C 

XT it > C7 — I !T0T!) C + 2| Uott> c 



\b 3 ) = ^(|omt> c -|otm>c+l itoit> c -| iron) 

l^> = ^(2| I TIT > c - I I ITT > a - I I TTI > c 



c 



-2|o|Tit> c +|o|in> c +|o|m> c ) 
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Id) 



^(|I010T> C + |0I0IT) C ) 



l<3»> = ^(|IloTO> c +|oolU> a ) 

Is.) = ^(inoor> c + |oonr> c 
\ c *) = ^(inoro> c -|oom> c 



where \abcde^ c refers to the cyclic state 



abcdej = —j=(\abcdej + \bcdea} + \cdeab} + \deabc) + \eabcd} 



Since the number of sites is odd, we can reproduce the BDS result with 
Applying (j2.1|) to (j2.13|) . we find the Hamiltonian matrix to be 



H 
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2g 
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-g 


-g 
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\l%9 
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-g 




















4 











-VSg 









-g 





-g 





3 


















-g 














3 





















-g 
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V2g 




V 














-V3g 








V2g 


3 


/ 



(2.13) 



(2.14) 

o ng. 



(2.15) 



Remarkably, the secular equation factorizes into a cubic and a sextic polynomial, with 

1 



det(A - H) 



P(A)Q(A) 



(2.16) 



where 



P(A) 
Q(A) 



A 3 - 11 A 2 - 10# 2 A + 39A + 34c/ 2 - 45 

4A 6 - 88A 5 - 30# 2 A 4 + 80A 4 + 432# 2 A 3 - 3904A 3 + 25# 4 A 2 
-2322# 2 A 2 + 10624A 2 - 170# 4 A + 5520# 2 A - 15360A 
+289# 4 - 4896/ + 9216. 



(2.17) 
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The desired solution is a zero of the cubic, which can easily be solved as a series expansion 
in g 2 , given by 



A 



5 + V - 6# 4 + 17 g* 



115# 8 833# 10 QU7g 12 44561c/ 14 303667# 



,16 



1778945^ 
64 



18 



2 4 
6255235# 20 
128 



16 



32 



45306735c/ 
256" 



22 



+ . . . . 



(2.18) 



This reproduces the result in ^3] using the Bethe ansatz. 

We can also write down a simple closed form expression for A in terms of A, where we 
find 



A 



57T 



(A 3 - 11A 2 + 39A - 45) 
(10A - 34) 



(2.19) 



From this we then deduce that for large A, A « V5A /(2n), and so using the dispersion 
relation in (J2.12)) . the up spins have momenta k = — 2tt/5, 0, 27r/5 and the down spins 
have k = —2ir/5, 2tt/5. Interestingly, for large A one sees that the A for this operator will 
be less than the A for the bare dimension 4 operator previously considered. 

The cubic -P(A) is itself the secular equation for a reduced Hamiltonian given by 



5 2y/2g 
H,, = I 2y/2g 3 y/2g 
V2g 3 



(2.20) 



while the sextic equation is the secular equation for 



( 








9 



\ 



3 9 4 Jig 



b 4 



o 






b 4 



-g 
-g 



o 





-g 



-g -g 

o o 

3 

3/ 



(2.21) 



The unitary transformation that rotates ()2.15|) into the block form H3 © Hq is a com- 
plicated g dependent transformation. This block diagonal form is probably due to the 
underlying integrability of the theory. 
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3 Large coupling for general states 



At large 't Hooft coupling A, the Hubbard model essentially reduces to a theory of free 
lattice fermions. Hence, the energies of the Hubbard states in the leading order approx- 
imation are explicitly known. What is not so trivial is determining which free fermion 
states go with the gauge invariant operators. In this section we will study this problem for 
the operators with two impurities and then provide some of the details for more general 
operators. Along the way, we will conjecture which Hubbard states correspond to gauge 
invariant operators for any value of A. 

We assume that L = 2K + 1 to avoid any phase 0. It is known from the 1 loop results 
that if M = 2 there are K possible solutions for the gauge operators, so we should expect 
the same number of solutions to apply in the free fermion limit of large A. One of these is 
the chiral primary which has A = L and has all momentum modes singly occupied. The 
operators which are not chiral primaries will have L — 2 different momentum modes for 
the up spins and two different modes for the down spins, with an overlap of the up and 
down modes. The sum of the momenta should be a multiple of 2n in order for the state 
to have cyclic symmetry. 

The problem now is to determine which momentum modes are occupied. To answer 
this, we consider the Bethe ansatz results in jT3]. The half-filled Hubbard model with 
= has a nested Bethe ansatz whose equations are [T7| 



n' 1 Uj - V2gsinq m -i/2 
Z—^f— ,,/o ' rn = 0...L 



j=1 Uj - V2g sin q m + i/2 



-|-r uj - vjty bin i; m tt/^ = tt uj ~ »fc -r o j = 1 M (3 1) 

*} Uj- V2g sin q m - i/2 ^1 u . _ Uk _ i ' 

The number of down spins is M and the energy of the state is given by 

L 

A = L + \^2g ^ cos q m . (3.2) 



m=l 



In the limit where g — > 0, the first Bethe equation implies 

e iq m L = e iq m ,L (33) 

for all m and m', while the second equation reduces to the Heisenberg chain Bethe equa- 
tions 

n^n;^f (3.4) 

m=l J kjtj J 
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The values of q m that correspond to the gauge invariant operators in the g — > limit have 
q m = 2nm/L for < m < L — 1 [TH], and so 

1 = TT«Lli/?. (3.5) 

Let us consider the case where M = 2. It was observed in ^Hj that the Bethe equations 
have a Z 2 symmetry under q m — > — g m , — > — In the Heisenberg limit, (|3.5|) implies 
M x = — m 2 = 11 > 0, so the corresponding solutions will maintain the Z 2 symmetry for 
g 7^ 0. Since L is odd, this means that go = for all g and q m = —qL- m - Hence the 
equations in (|3.1j) reduce to 



e^ L = " a -(v^™fln + i/2) a se ^ m = i...^ (3.6) 

u 2 — ( v 2(7 sin g m — i/2) 2 



TT ( M + V2) 2 -2^ 2 sin 2 g m = A i6m 
ll( M -V2) 2 -2^sin 2 g m -ll e " 



m=l 



At (7 = 0, we set the phases to Xm — and 8 m = 2im/K, for some n satisfying 1 < n < 
K/2. As g is turned on, u must start increasing, otherwise all 9 m will start decreasing 
and the equality in ()3.7j) will not be maintained. It is also clear that some 9 m must start 
decreasing toward while others will start increasing toward 2n, where those heading 
toward zero have the smaller values of sing m . In the limit of large g, it must happen that 
u ~ g, and that n of the 8 m approach 2tt. 

Examining (|3.6|) . one sees that Xm approaches if 6 m approaches 0, since the real 
parts of the numerator and demonitor in (|3.6|) differ from their counterparts in (J3.7|) by 
a finite amount, while the imaginary parts of the numerators and denominators do not 
change sign. This means that q m flows back to 2imi/L when g —>■ 00. However, Xm flows 
to — 2ir if 6 m flows to 2n, meaning that q m shifts to 2ir(m — 1)/L. If q m shifts and q m -i 
does not, then it must be that u flows to \/2g sin 27r (™~ 1 ) j n such a wa y that 

v 2 - 2g 2 sm 2 q m _ 1 > +1/4 
u 2 2g 2 sin 2 q m < 1/4 

In this analysis, we have assumed that q m is real throughout the flow. This is certainly 
true for g small enough. However, two q m can start moving into the complex plane once 
their real parts are equal. But as we have just seen, q m can equal q m -i only in the limit 
g —>■ 00, hence the q m remain real. 
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Figure 1: A flow with n = 2 for M = 2 and L — 17. After the flow there are doubly 
occupied modes at q = ±67r/17 and empty modes at q = ±107r/17. 



The largest values of sing m are centered around m = K/2. In the flow, there will be 
n sequential q m centered around m = K/2 that shift down to 27r(m — 1)/L. The negative 
q m of this sequence also shift. Hence, the final state for large g has a two doubly occupied 
modes at q = ±n(K — n + (3)/L and two unoccupied modes at q — ±tt(K + n + (3)/ L 
where (3 = (1 — (— l) K+n )/2. (Figure 1 shows the flow for L = 17 and n = 2.) Using 
perturbation theory about these free fermion states, the energies are found to be 

A - L + ^ (c OS p*Z?±>) - c OS (^±^») ) - ^ + O ( ±) . 

(3.8) 

Examining ()3.8|) , we see that L / vA is the effective expansion parameter and so (|3.8J) is 
a good approximation if L << v^A. In fact, the convergence appears better as n increases. 
In any case, L ~ V^X is where A crosses over from BMN to free lattice fermion behavior. 
In the appendix we compute the first order correction for all n as well as the second order 
correction for n « L. 

In principle this analysis can be done for M > 2. For one thing, it appears that the 
equations in (J3.1)) are consistent if all q m are real, as long as for every Uj there exists a 
Uk such that Uk = u*. The q m then remain real throughout the flows and no two q m will 
be equal for finite g, preventing any flow into the complex plane. In fact, what seems 
to distinguish between Hubbard states that correspond to gauge invariant operators and 
those that do not are that the former have real q m with ^2 Qm — mod 2n. To see that 
the latter do not have both properties, note that in the limit g — > 0, a state with total 
momentum zero but not corresponding to a gauge invariant operator will have A < L. 
The only way that this can occur is if cosg m ~ g^ 1 for some q m . But this itself only 
happens if Im q m ~ ± In g. 

The q m can be thought of as momentum modes for charge density waves, while the Uj 
are momentum modes for spin waves. Real q m means that no two charge density modes 
are bound together. In the g — > limit where the potential term dominates, this means 
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*H \ m*t H t H I H 



Figure 2: A flow with ni = 2 and ri2 = 1 for M + = M_ = 2 and L = 17. After the flow 
there are doubly occupied modes at q = ±67r/17 and q = ±8n/17 and empty modes at 
q = ±10tt/17 and q = ±12tt/17. 

that no two sites are doubly occupied. 

If the Uj are real, with M + of the positive and M_ negative, then our general 
expectation is that at least one of the positive q m shifts by —2itM + /L and at least one 
of the negative q m shifts by +2ttM^/L. It is also clear that if q m shifts by —2n£/L, then 
q m -\ and q m -2 must shift by at least — 2n(£ — 1)/L to avoid triple occupancy or modes 
hopping over each other. In general, for the positive q m we will find a set of sequences, 
where the first is a sequence of n\ modes shifting by —2tt/L, followed by a seqence of n 2 
shifting by —4ti/L, etc., all the way up to a sequence of hm + , where n« > 2 for i < M + . 
Behind this last sequence, there can be sequences shifting by ever smaller amounts. Figure 
2 shows the example with n\ = 2 and n 2 = 1 for M = 4 and L — 17. If some of the Uj 
are complex, then it will be possible to have multiple sets of these sequences. It is clear 
from this construction that the final configuration will have M of the momentum modes 
doubly occupied and M empty. 

As an example of this construction, consider the antiferromagnetic state (M = K) 
studied in [T51 [TH1 - This state would be expected to have all modes with < \k\ < tt/4 to 
be doubly occupied and those with 7r/4< \k\ <7r/2tobe empty. The sequences of shifts 
has rti = 2 for i < K/2. The total energy of this state at strong coupling is 

A AF w Ay/2g Y] cos ^ w L»l (3.9) 

n=l 

which agrees with the result in [T7j] . 

4 Comparing strong coupling to string theory 

Let us now ask if the Hubbard model at large A can reproduce nonclassical string results. 
At first sight this would seem problematic since generic string states are expected to have 
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energies satisfying 

A 2 -p 2 ^ An\f\ (4.1) 

for large A, where n is the level and p is a vector of Lorentz and i?-charges corresponding 
to the string center of mass motion This behavior can also be seen using a string 
Bethe ansatz [21]. But it seems from our discussion so far that A 2 will be proportional to 
A. In fact we can see that for any length chain and for any generic state of that chain, A 
will scale as vA for large enough A. As was emphasized in [T^], the Hubbard model does 
not have order of limits ambiguities, so any discrepancy cannot be due to this. 

In any event, let us see how one could possibly reproduce the string results by con- 
sidering the case M — 2, and assume that L » 1. In this limit (j3.8J) simplifies to 

A^L + — ^— - 2 (4.2) 

for n « L. The i?-charge of our state is J = L — 2 which corresponds to an angular 
momentum on S 5 . Hence we have 

A 2 - J 2 ss AnVX (4.3) 

which agrees with (|4.1jl if we identify n with the level n. Note that this approximation is 
valid as long as vA » L » A 1 / 4 , which is outside the BMN limit and thus beyond the 
range of validity for the BDS Bethe equations. Note that for a given n, if we reduce A 
such that we approach the BMN limit, then n is the level of the string in the plane wave 
limit [Hj. 

For states with M > 2, the energies will still have the form in (|3.8jl . that is, as a 
sum over doubly occupied modes minus a sum over empty modes. Hence, if the mode 
numbers are close enough to K/2 then A will still have the form in (|4.2|) for some integer 
n. However, it also seems that the lattice fermion levels do not precisely match with with 
the BMN levels. For example, for weak coupling and M = 4 we can build the BMN like 
state where U\ = — u 3 ~ L/(2tt) and u 2 = —u^ ~ L/(4tt). We would identify this state 
with a string state at level 3 0j. However, for strong coupling this state would flow to a 
state with two consecutive doubly occupied modes followed by two empty modes. Hence, 
this state has n = 4. It would be interesting to clarify this point. 

If L ~ A 1 / 4 then the approximation in (J4.3|) starts breaking down. If L « A 1/4 , then 
A 2 - J 2 » V\ for A >> 1. One possibility is that the string dual is a highly excited 
state, with excitations along all directions of AdSc, x S 5 . The Konishi operator would be 
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one example of this 2 . It is also possible that these states are dual to semiclassical long 
strings [T^l [201 [22] , which typically scale as y/X. 

There should be string states with small -R-charges and squared masses of the same 
order as that in (|4.1j) . but presumably their corresponding operators are not in the SU(2) 
sector. But it also seems that these operators might not have small bare dimensions 
either. Perhaps recent work on quantum integrability of the sigma model will lead to 
insights on this [23 [21 12S]- 

It is not clear how our results relate to the three loop discrepancy between gauge and 
classical string results j2El UH] • It was pointed out in [2Z| that the classical string solutions 
have small mixings outside of the SU(2) sector that could account for the discrepancy. 
This might require a generalization of the Hubbard model to include the entire operator 
spectrum in order to fully sort this out. One might make progress here by generalizing 
the antiferromagnet solution of ^3 fTU] . as was recently done for the SU(1, 1) sector in 
[2H [211 by using the asymptotic Bethe equations in HO] • 

Acknowledgments: This work was supported in part by Vetenskapradet. 

Appendix A: First and second order perturbative cor- 
rections 

We wish to find the first and second order perturbative corrections to states of the form 
|nj, — nj;n , —no), where ±nj are doubly occupied modes and ±n are empty modes. If 
m = n then the empty mode and the doubly occupied mode have combined to make 
singly occupied modes. Defining the Fourier transform of 




(A.l) 



the four fermion interaction becomes 



(A.2) 



Hence for the desired states, the first order correction to the energy is 

2(L — 2) 

A (i) = (nj, -nj,n , -n Q \H A \ni, -nj,n , -n ) = j • (A. 3) 



2 It has been been shown that in order for the Konishi operator to have A 1 / 4 behavior at strong 
coupling, it cannot be a rational function of y/X |2()| . 
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For the second order correction, we assume that nj = (K — n + /3)/2 and n$ = 
(K + n + /3)/2, where f3 is defined in the text. In general, we have the following matrix 
elements 



(nj, — nj,n Q , — no | #4 1 raj + m, -nj; n + m, -n Q ) 
(nj, — nj, no, -n |-ff 4 |nj, — nj + m; n , -n + m) 
(nj, -nj,n , -n |i/" 4 |^j + m , n , -n + m) 
(nj, — nj, n , — no|# 4 pj, — + m;n + m, -n ) 

We also have that the difference in energies are 
A (o) _ A (o) 

n-t ,— nj ,no,— no n|+m,— nj ;no+m,— no 



-1/L 
-1/L 
-1/L 
-1/L. 



y/X 
2vr 



n(K-n + /3) n(K + n + {3) 

cos — cos ■ 



— COS 



2K + 1 

ir(K - n + 2m + 0) 
2K+1 



2K+1 

n(K - n + 2m + (3) 



cos 



2K + 1 



rnr 2 y/\ 



L 3 



(2m 2 + (1 - 2/3)m) 



A (o) _ A (o) 

,— n-j ,no,— no n|+m,— nj;no,— no+m 



2vr 



cos 



2K + 1 



cos ■ 



ir(K + n + (3) 
2K + 1 



7r(K - n + 2m + (3) ir(K - n + 2m + (3) 

COS t-tz : h COS — 



2K+1 



2K+1 



(2m) 



(A.4) 



(A.5) 



(A.6) 



Since (|A.5|) is large compared to (|A.6|) . we can ignore the contributions coming from the 
nj + m, — nj; n , — n + m) and |nj, — nj + m; n + m, — n ) states. We then find, 

L 3 (-l/L) 2 4L 



A (2) 



nvr 2 v / A(2m 2 + (1 - 2/3)m) nvr 2 ^ 



(A.7) 
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